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Partitioned block frequency domain adaptive filter 



The invention relates to a partitioned block frequency domain adaptive filter 
for filtering an input signal in dependence on a control signal, the adaptive filter comprising a 
plurality of parallel arranged filter partitions, each filter partition being arranged for modeling 
a part of an impulse response of the adaptive filter, each filter partition having update means 
5 for updating filter coefficients of that filter partition by circular convolutjng a signal 

representative of the input signal and a signal representative of the control signal, the update 
means comprising constraint means for intermittently constraining the filter coefficients by 
eliminating circular wrap-around artifacts of the circular convolution. 

The invention further relates to an acoustic echo canceller comprising such a 
10 partitioned block frequency domain adaptive filter and to a method of adaptively filtering an 
input signal in dependence on a control signal. 



A partitioned block frequency domain adaptive filter (PBFDAF) according to 
15 the preamble is known from the paper "Multidelay block frequency domain adaptive filter", 
IEEE Transactions on Acoustics, Speech and Signal Processing, vol. 38, no. 2, pp. 373-376, 
by J.S, Soo and I£.K. Pang. Adaptive filters can be very useful devices in several applications 
of digital signal processing. Examples of such applications are: channel equalization, array 
signal processing and noise and echo cancellation. Here, the Acoustic Echo Canceller (AEC) 
20 is merely used as a means for demonstrating the invention, and it is noted that the invention is 
not restricted to this specific application. Acoustic Echo Cancellers are used, for example, in 
teleconferencing systems or in speech recognition systems. The structure of an AEC is 
depicted in Figure I and will be explained as part of a teleconferencing system. 

In a teleconferencing system with two speakers in two different rooms, an . 
25 AEC is implemented twice. The speech signal from the far-end speaker x[k] is generated in 
the near-end room by the loudspeaker 14. The transmission between the loudspeaker 14 and 
the microphone 16 can be modeled by a room impulse response. The echo signal e[k], picked 
up by the microphone 16, can be seen as a convolution of the signal x[k] with this room 
impulse response. In the case that both parties are speaking (i.e. double-talk), the microphone 
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16 also picks up a near-end speech signal s[k]. If there was no adaptive filter present, the 
signal 

would be directly passed on to the far-end speaker. Besides the wanted near-end speech 
signal, the far-end speaker also hears echoes of his own speech signal. For the overall 
5 teleconferencing system, this results in a non-optimal communication between the two 
speakers. 

The adaptive filter, which is formed by the update means or correlation means 
10 and the convolution means 12, reduces these unwanted echoes by modeling the real room 
impulse response with a Finite Impulse Response (FIR) filter. The filter-coefficients or filter- 

10 weights w[k] of this FIR filter are updated in the update means 10 by correlating the residual 
signal r[k] with the input signal x[k]. By convoluting or convolving the input signal x[k] with 
the filter-coefficients w[k] in the convolution means 14, the adaptive filter estimates the 
unknown acoustic echo signal e[k], indicated by £[k]. This estimate of the echo is subtracted 
from the real echo by means of a subtracter 18, As a result, the echo level is reduced, leading 

IS to improved communication between the speakers. 

For efficient implementations of adaptive filters, it is known to use block 
signal processing in combination with Fast Fourier Transforms (FFTs) for performing fast 
convolution (filtering) which permits adaptation of filter parameters in the frequency domain 
in a computationally efficient manner. To do this, a block of input-samples is collected and 

20 the adaptive filtering is performed in frequency-domain. Commonly, Fast Fourier Transforms 
(FFTs) are used to calculate the frequency- domain data from the time-domain data although 
it is noted that also other transforms (e.g. (I)DCT transformations) can be used for this 
purpose. 

An elementwise multiplication of two frequency-domain vectors corresponds 
25 to a convolution in the time-domain. If both vectors are transformed with the Fast Fourier 

Transform this convolution is a circular convolution. In adaptive filtering linear convolutions 
in stead of circular convolutions are needed to convolve an infinite length input data stream 
with a finite length vector. By means of the so-called Overlap-Save (OLS) method a linear 
convolution can be performed using FFT/IFFTs. It is noted that the Overlap-Add (OLA) 
30 method can also be used for efficiently convolving an infinite length vector with a finite 
length vector. However, adaptive filters using the OLS are usually less computational 
complex than adaptive filters that use the OLA method- 
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An adaptive filter comprises a linear correlation part 10 and a linear 
convolution part 12. The convolution part 12 is needed for filtering an input signal with the 
filter-weights, while the correlation part 10 calculates a gradient estimation needed to update 
these adaptive filter-weights. After convergence the adaptive weights have adjusted to the 
5 optimum solution and represent the room impulse response. The correlation part 10 and the 
convolution part 12 may both be implemented by means of the Overlap-Save method. The 
correlation part 10 is similar to the convolution part 12 in that in both parts a linear 
convolution is calculated. The correlation part 10 differs from the convolution part 12 in that 
in the correlation part 10 one input block has to be mirrored in time. This mirroring in time . 
1 0 can be performed in frequency domain by means of a complex conjugate operator 32 as 
shown in Figure 2. 

A drawback of (50% overlap) block frequency domain adaptive filters 
(BFDAFs) is the delay of N samples due to the block processing. First a complete block with 
input samples is collected, before the block can be processed. This problem is solved in 
1 5 PBFD AFs where this block size is reduced by partitioning of the adaptive filter. A PBFDAF 
can be seen as a parallel implementation of a number of BFDAFs which all model a small 
part of the impulse response. The convolution result is obtained by adding the results of all 
the convolutions. 

Figure 2 shows a block diagram of a PBFDAF. Samples of the input signal 
20 x(k] are collected into blocks with 50% overlap by means pf the serial to parallel converter 
20 and the 50% overlap block 22, These blocks of input samples are then converted into the 
frequency domain by a FFT 26. In block 42 N zeros are appended on the left of a block of N 
samples of the residual signal or control signal r[k]. The resulting block of 2N time domain 
values is thereafter converted into the frequency domain by another FFT 26. Next, the power 
25 of this frequency domain signal is normalized by means of an elementwise multiplication 
with a normalization factor in an elementwise multiplier 40. 

Each partition comprises a complex conjugate operator 32, an update block 34 
and an elementwise multiplier 40. Additionally, all partitions with the exception of the left- 
most partition, cgmprise a block delay element 30. The update means or correlation means 10 
30 comprise the complex conjugate operator 32 and the update block 34, The convolution means 
12 comprise the elementwise multiplier 40. The output of the complex conjugate operator 32 
is the signal representative of the input signal. The output of the elementwise multiplier 40 is 
the signal representative of the control signal. The update block 34 updates and outputs the 
filter coefficients on the basis of the signal representative of the input signal and the signal 
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representative of the control signal. Next, the updated filter coefficients are convoluted with 
blocks of samples of the input signal by means of the elementwise. multiplier 36. The outputs 
of the elementwise multipliers 36 of all the partitions, i.e. the results of all the partial 
convolutions, are added together in an adder 38 in order to form the overall convolution 
5 result. This overall convolution result is then transformed into the time domain by means of 
an IFFT 28. FinaHy, the left-most N samples of the block of 2N time domain samples of the 
overall convolution result are discarded in order to form the estimated echo signal S|k]. This 
estimated echo signal is subtracted from the signal picked up at the near end by means of 
subtracter 1 8. The result of this subtraction is the residual signal r[k]. 
1 0 The update means comprise the complex conjugate operator 32, the update 

block 34 and the elementwise multiplier 36. 

Figure 3 shows a block diagram of an embodiment of the update block 34. m 
this embodiment the adaptive weights are stored in time-domain. Note t ha t the index i 
indicates the partition number. Furthermore, the following variable W* N is defined as; 



-^(^ v ) 4*w?IkN}- (3.1) 



The signal representative of the input signal and the signal representative of the control 
signal are elementwise multiplied with each other by elementwise multiplier 50. This 
elementwise multiplication in frequency domain constitutes a circular convolution in time 

20 domain. The result of this convolution comprises circular wrap-around artifacts which are 
removed by means of a rectangular time domain constraint window which is formed by IFFT 
52, element disposal window 54, adder 56, block delay element 58, zero appender 60 and 
FFT 62. First, the convolution result is transformed into the time-domain by the IFFT 52. 
Next, the window 54 discards N points of the convolution result. This is needed to remove 

25 the circular wrap-around artifacts from the Overlap-Save method. The remaining N samples 
are used to update the adaptive weights with the adder 56 and the delay element 58. After this 
update, N zeros are augmented. This is needed for the efficient convolution of the adaptive 
weights with the input signal, by means of the Overlap-Save method. Finally, the result is 
transformed back Jp frequency-domain. 

30 The two mirroring operators J N can be eliminated by storing the vector 
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in time-reverse order. Furthermore, the two time-domain windows, the FFT and the IFFT, 
can be combined in a tune-domain gradient constraint operator. As a consequence of the 
gradient constraint operator or constraint operator, the adaptive coefficients need to be stored 
_ and_up-dated_ in frequencyjiomain instead of time-domain. These frequency-domain 
5 coefficients will be denoted by N . The separation of the constraint operator and the 
adaptive coefficients can be illustrated by the following deduction: 

H$»i<* + DAT) » ? ( ^ ) J* (i^VATl + J* ( I* O v ) F^Z^lxN]) 

« fgHKNl+J^^^Z^lKNl* (3.2) 

where 

10 is the rectangular constraint window vector, containing ones in die first N elements and zeros 
in the last N elements. Here, this constraint window is referred to as the rectangular 
constraint window. The constraint means comprise this rectangular constraint window. 
Figure 4 shows the frequency-domain coefficient update together with the rectangular 
constraint window. 

IS It is to be noted that there are reasons for updating and storing the adaptive 

weight in frequency-domain, despite the fact that frequency-domain weights need more 
storage elements compared with the time-domain weights. The most important reason being 
the possibility to omit the gradient constraint, in order to reduce the computational 
complexity. 

20 Mostly, the rectangular constraint is placed before the update, but it is 

advisable to place this constraint after the update as is shown in Figure 5. Calculating the 
rectangular constraint with infinite precision gives results which are the same for both 
situations (Figures 4 and 5). However when calculating the weight update in Figure 4 with 
finite precision (e.g. fixed-point arithmetic), numerical round-off errors are accumulated in 

25 the weight update. Eventually, after many iterations of the weight update, this leads to large 
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6 21.08.2000 
errors in the last N points of the coefficients- Correction of these accumulated errors is 
prevented by the right half part of the constraint When placing the constraint after the weight 
update as depicted in Figure 5, the result after calculating the constraint also contains 
numerical round-off errors in cases of finite precision. However, as the constraint operation is 
5 performed on the actual coefficients, there is no accumulation of these errors which could 
lead to long-term problems. Hence, rectangular constraint operations are preferably placed 
after the weight update. 

Figure 6 shows a block diagram of an embodiment of the update block 34 
which illustrates the operation of the PBFDAF which is disclosed in the above mentioned 
10 paper. The use of the rectangular constraint window 82 is computationally relatively 
complex. This complexity can be lowered by constraining the coefficient vector of a 
particular partition only once in a while. For this purpose the update block 34 comprises two 
switches 80 and 84. The operation of these switches 80 and 84 may be described by the 
following function: 



f ) i 
\ 0 c 



if k mod KP »i 

Otherwise ' t3 **> 
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In this formula the parameter K indicates the number of partitions and the parameter P is the 
alternative constraining period in number of block iterations of the algorithm. The 
rectangular constraint window is calculated only when the switches 80 and 84 are both in the 
20 1 -position. When die switches 80 and 84 are both in the 0-position the convolution results are 
not constrained at all. 



It is an object of the invention to provide a PBFDAF with an improved 
25 convergence behavior, while maintaining a relatively low computational complexity. This 
object is achieved in the PBFDAF according to the invention, which is characterized in that 
the update means further comprise selection means for selecting and removing at least part of 
the circular wrap-around artifacts, the selection means comprising an approximation of a 
rectangular constraint window. Because the constraining of the circular wrap-around artifacts 
30 is performed intermittently, i.e. only once in a while, it is possible that a too large 

accumulation of circular wrap-around artifacts takes place in between these constrainings, 
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7 21.08.2000 
resulting in a relatively bad convergence behavior. This is prevented by applying an 
approximation of a rectangular constraint window for the purpose of selecting and removing 
at least part of the circular wrap-around artifacts. The approximation of the rectangular 
constraint window allows to eliminate a large part of the circular wrap-around artifacts in a 
5 _ computationaUy.efficient way. 

An embodiment of the PBFDAF according to the invention is characterized in 
that the selectiomneans are arranged for selecting the at least part of the circular wrap- 
around artifacts when the filter coefficients are not constrained by the constraint means. By 
this measure the selection and removal of the at least part of the circular wrap-around 

10 artifacts is not performed when it has no use, i.e. on the moments that the circular wrap- 
around artifacts are eliminated anyway by the constraint means. 

An embodiment of the PBFDAF according to the invention is characterized in 
that time domain values of the approximation are larger than or equal to zero. Negative 
values in the time-domain window need to be avoided in order to avoid divergence of the 

15 adaptive weights. The reason for the divergence in the case of negative values in the 

constraint window can be explained in an intuitive way, Because of the negative values, the 
correlation part produces results that are negatively correlated with the residual signal. Hence 
the adaptive weights are updated in die wrong direction. As the adaptive weights become less 
accurate, there is more negative correlation and worsening of the weights in every next 

20 block-iteration. 

An embodiment of the PBFDAF according to the invention is characterized in 
that the approximation in time domain has substantially high slopes near the positions which 
correspond to the positions of the transitions in a rectangular constraint window. By this 
measure, a relatively large amount of the circular results that converge to a non-zero value is 

25 eliminated by the selection means. For example, when using a 50% overlap, the high slopes 
Should be positioned at 0 and N, with N being the block size of the filter and 2N being the 
FFT size. When using a 75% overlap, the high slopes should be positioned at 0 and 3N, with 
N being the block size of the filter and 4N being the FFT size. 

An embodiment of the PBFDAF according to the invention is characterized in 

30 that frequency domain values of the approximation each comprise a real value and complex 
conjugate imaginary values, whereby at least part of the imaginary values form a row of 
numbers, the numbers being obtainable from one another by multiplication. By this measure, 
at least part of the frequency domain values of the approximation can easily be calculated, i.e. 
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8 21.08.2000 
merely by multiplying other frequency domain values with a multiplication factor, resulting 
in a PBFDAF with a relatively low computational complexity. 

An embodiment of the PBFDAF according to the invention is characterized in 
that the frequency domain values of the approximation are defined as: 



2AT<? for i=0 

0 for 0 < ; <2^i even 



- «?"(^ 1=t )^ j for 0 < i < 2N, i odd 

with / being an index number, m being a multiplication factor, a being a mean value. Because 
of the relation between two subsequent frequency components in this constraint-window, the 
frequency domain circular-convolution of this window with a complex valued Hermitian 
signal can be calculated in a very efficient way. 



10 



The above object and features of the present invention will be more apparent 
from the following description of the preferred embodiments with reference to the drawings, 
15 wherein: 

Figure 1 shows a block diagram of an embodiment of the acoustic echo 
canceller, w 

Figure 2 shows a block diagram of an embodiment of a PBFDAF comprising 
an update block 34, 

20 Figures 3 to 7 show block diagrams of embodiments of the update block 34, 

Figures 8 to 13 show diagrams by means of which the compensation of 
circular wrap-around artifacts in the PBFDAF according to the invention is illustrated, 

Figure 14 shows a block diagram of an embodiment of a generator means for 
generating forward and backward compensation signals, 
25 Figure 15 shows some signal diagrams by means of which the operation of the 

generator means as shown in Figure 14 is explained, 

Figure 16 shows a block diagram of another embodiment of a generator means 
for generating forward and backward compensation signals, 

Figure 17 shows some signal diagram s by means of which the operation of the 
30 generator means a& shown in Figure 16 is explained, 
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9 21.08.2000 
Figures 1 8 and 1 9 show some signal graphs illustrating the operation of the 
generator means, 

Figure 20 shows a block diagram of an embodiment of a PBPDAF according 
to the invention, 

5 Figure 21 shows a block diagram of an embodiment of the update block 34 for 

use in the PBFDAF of Figure 20, 

Figure 22 shows a block diagram of a further embodiment of a generator 
means for generating forward and backward compensation signals, 

Figure 23 shows a block diagram of an embodiment of a generator means for 
1 0 generating a backward compensation signal. 

In the Figures, identical parts are provided with the same reference numbers. 
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15 In the Figures and in this description uppercase symbols are used to denote 

frequency-domain variables and lowercase symbols for time-domain variables. A boldface 
font is used to denote matrices, while vectors are denoted by underlining fonts. Furthermore 
the superscripts t , * and h are used to denote the transpose, complex conjugate and the 
complex conjugate transpose (Hennitian) respectively. O 2 * and I N represent a zero and 

20 identity matrix, the superscript indicating the dimension of the matrix. Generally, superscripts 
indicate the dimension and subscripts indicate the element number of the vector or matrix. As 
most of the vector^ and matrices used here have dimension 2N x 1 and 2N x 2N, a vector or 
matrix with omitted superscript indicates a vector of length 2N or a matrix with dimension 
2N x 2N- The indices between square brackets (e.g. x[k]) s indicate a time-index, diag {.} 

25 denotes the operator that transforms a 2N x 1 vector into a 2N x 2N matrix containing the 
vector coordinates as diagonal elements while other elements are zero. The mathematical 

expectation is denoted by 

A circular shift of the data over L positions in an 2N dimensional vector is 
carried out by the matrix; 

30 
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which the zero matrices O have appropriate dimensions. Note that 

The data of an 2N dimensional vector can be mirrored with the 2N x 2N mirrored matrix J 2 * 
that is defined as; 



A reverse circular shift over L positions, in opposite direction to 
is carried out by the matrix 

Note that 
J§* J 2 * 

1 0 Besides mathematical descriptions also figures based on signal processing 

blocks are used. The explanation of these blocks is appended to the description. There also an 
overview is given of the symbols used. It is to be noted that although in this description only 
50% overlap or Half Overlap (HOL) of block-data is used, the invention is not limited to this 
specific implementation of overlap. The invention can also be used with other overlap 

15 percentages. 

The Normalized Mean Square Error (NMSE) 3 in dB, is a performance measure 
of adaptive filters and is defined as follows: 
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MMSEfAj » 10 log j . (1.7) 

where N is the length of the block that contains the residual samples for which we calculate 
the NMSEfk]. The steady-state NMSE is the value of the NMSE after convergence of the 
adaptive filter, and is a measure of the residual echo. The lower the steady-state NMSE, the 
better the performance of the adaptive filter. 
5 The CTC, which is another performance measure of adaptive filters, is defined 

as the slope of the NMSE curve, during the transient of the algorithm. It assumes that the 
initial part of the learning curve is linear- The CTC is defined in dB/ms. 

An EPFT and FFT operation and an elementwise multiplication are needed to 
calculate the constraint in order to null half of the coefficients in time-domain. This requires 

10 operations, where 2N is the length of the FFT/IFFT. This constraint can also be performed by 
calculating a circular convolution directly in frequency-domain. The straightforward 
calculation of the circular convolution requires 

operations. The frequency-domain coefficients 

for the circular convolution with the rectangular constraint are: 

15 

12N*\ fori-Q 
0 for 0 < i < 2N, i even (3 g* 

v ft>rO<f <2*.todd 

When using a first order approximation, only use the real coefficient and the 
first two mutually complex conjugates are used. This results in a sinusoid constraint window 
that can be computed in frequency-domain with much less computational complexity than by 
means of an FFT/IFFT, Off course, also other low-order approximations of the rectangular . 
20 constraint window could be used 

Constraint windows with other than 50% overlap (half overlap) contain 
frequency components that have both real and complex values. Hence such a calculation 
needs additional computational complexity compared with the 50% overlap windows. As 
mentioned before, in most practical cases half overlap constraint windows are used. 
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12 21.08.2000 
The rectangular constraints need to be placed after the update. This is however 
different for the constraint approximations. If the approximations were placed after the 
update, the time-domain adaptive weights would constantly be multiplied by the weights of 
the approximation window. As most of the weights in the approximation window are 
5 generally not equal to one, this would lead to adjustments of the weights in every block- 
iteration and the adaptive filter would not converge. Hence all constraint approximations are 
best applied before the update, e.g. before the adder 56. 

Negative values in the time-domain window need to be avoided in order to 
avoid divergence of the adaptive weights. The reason for the divergence in the case of 
1 0 negative values in the constraint window can be explained in an intuitive way. Because of the 
negative values, the correlation part produces results that are negatively correlated with die 
residual signal. Hence the adaptive weights are updated in the wrong direction. As the 
adaptive weights become less accurate, there is more negative correlation and worsening of 
the weights in every next block-iteration. 
15 As the frequency-domain coefficients of equation 3 . 8 result in time-domain 

coefficients that can be negative, we need to scale the frequency coefficients. Using more 
frequency components in the constraint approximation, results in better performance, but 
enlarges the computational complexity. To achieve an optimal performance with low 
computational complexity, it is investigated which properties of the constraint are important. 
20 Therefore we perform an experiment with two constraint windows. 

Constraint window 1 includes the real and first 3 mutually complex conjugate 
components according to formula 3.8. A scaling by a factor 5/6 is applied in order to avoid 
negative values in the window. The second window contains the same gradient in the slope 
compared with the first window type, but is more accurate (with respect to the rectangular 
25 constraint window) in between the slopes. The step-size parameter 2a is set to the value of 
0.5/K. Both constraint approximation windows lead to a less accurate steady-state NMSE 
compared to the full-constrained case. The difference between the two approximation 
windows however do not differ very much. Apparently, the values in between the slopes are 
of little importance for the convergence behavior. This can be explained due to the fact that 
30 the time-domain adaptive weights in between the slopes are updated with random (noisy) 
values. Hence the weights do not converge and are zero in average. This is not true for the 
circular artifacts closer to the slope. Hie adaptive weights close to these slopes converge to a 
non-zero value and introduce errors in the convolution part of the PBFDAF algorithm, as the 
neighbor partitions are coupled. To eliminate as much of the circular results that converge to 
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a non-2ero value, high slopes should be created on the two places where there are transitions 
in the rectangular constraint window. 

As high slopes are important for obtaining a good convergence behavior, as 
many frequency components as possible should be used. When using a regular circular 

5 convolution_algonthm 5 mcoiporatmg.more frequency components leads to a linear increase in _ 

computational complexity. Therefore a regular circular convolution algorithm should not be 
used. 

An efficient method for convolving a signal with an approximation of the 
rectangular con-straint is by assuming that the odd frequency components in equation 3.8 are 
1 0 related with the neighbor odd components by a multiplicative factor. As this window also 
contains higher-order frequency components, the slope of this constraint approximation is , 
very steep. Hence this window is called the high-slope constraint window. The high-slope 
constraint window CP 1 in frequency domain for 50% overlap partitions is defined as: 



(2?*), 



2N-a i =0 

0 for 0 < i < 2N r i even ^ 4 ^ 

_3^£ (^) - »r(*^)) j for 0 < i < 2U, i odd 

1 5 Because of the relation between two subsequent frequency components of this above defined 
constraint-window, the frequency domain circular-convolution of this window with a 
complex valued Hermitian signal can be calculated in a very efficient way. 

Experiments have shown that good results can be obtained when the 
constraint-window g 21 * in time domain qualifies the following conditions: 

20 

2. High-slppe on / = Oand/ » 'V, 

From the experiments it follows that a steep transition on i = 0 and 
i = N (i.e. high-slope) is of more importance than the nulling in the range N £ i < 2N 
By choosing m- 2.166 and a = 0,57, the high-slope constraint window 

just qualifies the first condition above, while maximizing the derivative of the slope in the 
25 middle, as proposed in condition 2. This optimal value of m is found experimentally. With 
the value a, the area between the zero level and the high-slope window can be adjusted. For 
a=0.57, the mean amplitude in the range 0 ^ i < N is normalized to 1 . 
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To derive an efficient algorithm for convolving the high-slope window with 
the complex-valued Hennitian input-signal, we start with the definition of the circular- 
convolution: 

■ . 2AT-I 

_ ®G>») k » J- £ ^) (t _,, Marf (g*) r (4.3) 

j"=0 

As the circular convolution is a linear operation, we can apply the superposition theorem and 
split up the circular convolution with G 2 * into three parts: 

(4-4) 



where: 

£/v #«M> f f add v z / 

* ,r f=0, J odd ^ * / 



As there is a direct relation between formula 4,5 and formula 4.6, we only need to cal- 
culate one of these circular convolutions in the algorithm. This relation between the two 
formulas is expressed in the following lemma: 

Lemma 4^1 if (G^) fc » (fif')L/-* ™* t* 2 *)* « U?*)w-* 
(?.<?. A; 2 ^ Hermtii<w), then: 

(X™ * G^) k - QP» ® flj*)^ (4.7) 



Hence, in the next derivation of the algorithm, we only consider the circular convolution of 
formula 4,5. Because each next convolution-element 

can be expressed as a function of the actual convolution-element 

a very efficient convolution algorithm can be derived. The iterative expression is formulated 
15 in the following lemma: 
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Lemma 4^2 if £f* ^ | " 

U^h - AT 2 " HermiOaa). Thar: 

^(^■(^•^.-<^W|-(-.----»-«)/). (4.8) 

On a similar way as the above derivation, we can derive an expression for calculating a 
convolution element 

as a function of the convolution element 

<s$*).- 

This results in the following expression: 

w "'(s7(^«^*')» + <^)*-.i('»-'»" w - ,> )-')- <"•»> 

5 Note that the above two expressions are numerically stable iterations, because of the 

multiplication with the factor m"\ One can also derive an expression where in each iteration a 
multiplication with factor m is done; however, such an expression leads to less accurate 
results. 

In order to calculate a convolution with the high-slope constraint in frequency- 
10 domain, we basically need to do the following sub-calculations: 

1. Calculation of ^ (X 2JV © a) (i.e scale operation). 

3. Iterative calculation of the convolution by means of lemma 4,2.2. 

4. Calculation of j^- @ £r*) k out-of the calculations above 

The convolution with the mean-value a can be calculated by scaling the input-signal X 2 * with 
the factor a. We note that when performing more than one circular convolution with the high- 
slope constraint window (as is the case in the PBFDAF algorithm), this convolution with the 
mean-value 

1 5 can be performed at the cost of only 1 multiplication. 
The calculation of 
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• ^ «^ 4r • fif*), 
is needed for the first iteration of the iterative calculation of the circular convolution, one for 
the odd points and one for the even points. As the sum of the first few mutually complex 
conjugate frequency components of equation 3.8 is almost equivalent to the total sum of all 
mutually complex conjugate frequency components, there is no need to calculate 

5 and 

based on all the coefficients of 

Therefore we perform the initial calculation based upon the following calculation: 

m^®s&\*±±&™), k ^ m mx <4.io) 

In the calculation of 

£<(^*®^^ 
1 or 2 coefficients of 

10 need to be considered (i.e. Q — 1,3) in order to get accurate results. 
The fector 

in lemma 4.2.2 can be precalculated beforehand by the input-data X 2 * . As the multiplication 
with a can be performed elsewhere, as described in 4.2.2, the precalculation factor is 
J. (« x 

For large N, the pre-calculation is approximately m/2. m ~ 2.166 results in a pre-calculation 
15 of 1 .083 . We can approximate this value by 1 A resulting in a reduced complexity of 2N 
multiplications, without violating the condition 

a. v. 

The complexity of the iteration-steps can be reduced even further, as a 
multiplication by two can be performed by a shift operation. When using a precalculation of 
20 1 .0 and m =2, the constraint-window slightly drops below zero, thus violating the condition 
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As negative values in the constraint should be avoided, it is not advisable to use m = 2, When 
we perform the iterative calculations with the pre-multiplication of 1.0, a single iteration in 
the convolution algorithm is done at the cost of two additions and one multiplication. 

Because of the relation between 
GfandCf 



5 according to lemma 4.2.1, only the convolution with G, need to be calculated. Also half of 
the output-points of 

need to be calculated, as the frequency-domain vectors are Hermitian. 

Figure 7 shows an embodiment of the update block 34 in which an 
approximation of a rectangular constraint window 92, e.g. the high slope constraint window, 

10 is employed in combination with the rectangular constraint window or constraint means 82. 
The selection means comprise in this embodiment the approximation of a rectangular 
constraint window 92. The switches 80 and 84 control the operation of the update block 34. 
When both switches 80 and 84 are in the 1 -position, the exact rectangular constraint window 
82 is used to eliminate all the accumulated circular wrap-around artifacts. This exact 

1 5 rectangular constraint has a relatively high computational complexity and is performed only 
once in a while. When both switches 80 and 84 are in the 0-position, the approximation of the 
rectangular constraint window 92 is applied in order to eliminate a large part of the circular 
wrap-around artifacts in a computationally efficient way. In this way it is prevented that a too 
large accumulation of the circular wrap-around artifacts takes place, which would result in a 

20 relatively bad convergence behavior. 

Experiments have shown that the PBFDAF as described in die above 
mentioned paper exhibits jumps in the convergence curve at the moments when the 
rectangular constraint 82 is applied: at these moments the NMSE is suddenly worse than 
before. Since the degradation in the NMSE occurs exactly after the moment that the 

25 rectangular constraint is applied, the removal of the circular-errors are responsible for this 
phenomenon. Apparently the accumulated circular-errors in the adaptive weights somehow 
contribute to the model of the room impulse response, and removing these circular-errors by 
means of rectangular constraining results in a temporarily degradation of the NMSE in the 
next block-iteration- 

30 We start with an analysis, where we demonstrate that it is useful to 

compensate some of the wrap-around artifacts. This analysis is done by means of graphical 
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18 .21.08.2000 
construction of the convolution part in the adaptive filter. After this analysis we show how 
the compensation signals can be constructed. 

We can eliminate the temporarily degradation in the convergence curve of the 
known PBFDAP by compensating the removed circular-errors in the the other partitions after 
5 a rectangular constraint is applied. As a result, the output of the partitioned convolution is 
approximately the same as when the constraint was not applied For explaining this method, 
we consider Figure 8 which shows the input-blocks in time-domain. After the circular 
correlation of the residual signal with the input signal blocks is performed, the result 

contains the wanted linear correlation in the range [0..N-1], while the range [N..2N-1] is 
1 0 polluted with circular wrap-around artifacts. In the full constrained PBFDAF all these wrap- 
around artifacts are eliminated. In this analysis we assume that the results near the correlation 
outputs N and 2N can be considered to approximate a linear correlation. In the PBFDAF 
without power normalization, this can be easily understood by correlating the residual signal 
in time-domain (augmented with N 2eros) with the input signal blocks of Figure 8. For the 
1 5 correlation output points higher but close to N, the amount of wrap-around artifacts is low. 
This is also true for the points lower but close to 2N- For correlation output points in the 
direction of 3N/2 the approximation gets worse as the amount of wrap-around artifacts 
increase. Therefore we refer the correlation output points higher but close to N as being the 
left tail and the points lower but close to 2N as being the right tail. As the right tail 
20 approximates the linear correlation for negative time lags, it is better to change the principal 
interval as is depicted in Figure 9. As both tails are approximately results of a linear 
correlation, the coefficients of the tails in the adaptive weights 

(in time-domain) converge to a value that is approximately correct. The tails of these 
adaptive weights can be successfully used for eliminating the temporarily degradation in the 

25 convergence curve of the known PBFDAF. On the other hand, the correlation results in the 
neighborhood of 3N/2 (with a high amount of wrap-around artifacts) are useless, due to the 
fact that there is an approximately 50% wrap-around between the two input blocks of length 
N. Hence, these correlation results have nothing to do with linear correlations. As these 
correlation results can be assumed random, the adaptive weights in the adaptive filter are also 

30 updated randomly, As a result, these parts in the adaptive weights are fluctuating and do not 
converge to a stable value, Therefore we refer these values in-between the two tails as the 
noisy part. 
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The circular wrap-around artifacts introduce errors in the convolution part of 
the adaptive filter, as neighbor partitions in this convolution are coupled- Because the two 
tails contain values that are approximate results of linear correlations, we can reuse these two 
tails on the moments where the rectangular constraint is applied in the known PBFDAF. This 
-5 - is done in such a way that the output of the block convolution in the adaptive filter remains ~ 
approximately unchanged, while all circular wrap-around artifacts (i.e, the two tails and the 
noisy part) in thatspecific partition are eliminated. A method for achieving this, is by 
compensating the left tail into the next partition and compensating the right tail into the 
previous partition. 

0 We explain the method of compensation graphically. In Figure 10 we see the 

two weight vectors of the adaptive filter. For the adaptive weight vector 

the left part represents the weights based on the correct linear correlations. The right side of 
this vector contains weights polluted with circular wrap-around artifacts. The weights are 
constructed out of the two tails as discussed earlier. The left tail is indicated with the number 
5 1 3 while the right tail is indicated with number 2. In between these tails, the weights are noisy 
and therefore useless. 

In the previous section we discussed the weights in the adaptive filter as being 
constructed out of results related with linear correlations and the circular wrap-around 
artifacts. It is also explained that a part of these wrap-around artifacts (the two tails) are 
useful results. The compensation discussed in this section explains how these tails can be 
reused. As we have two useful tails we can use them for a backward and a forward 
compensation. Figures 1 1 and 12 are used for explaining the forward and backward 
compensation. 

The circular convolution of 
xj#[K X} with v^{kN] 

can be performed graphically by mirroring one signal (e.g. the adaptive weights) and moving 
this signal to the right; this is indicated in Figure 1 1 . The total sum of products of both 
signals in one convolution-iteration is the output-element of that specific iteration. As we 
perform a circular convolution, shifting the adaptive weights is actually a rotation of the 2N 
length vector. The convolution output can be described by the following formula: 
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In Figure 1 1, three stages of graphical construction of the convolution-output are shown, We 
note that this figure shows the circular convolution in time-domain while in our efficient 
implementation, the circular-convolution is performed in frequency-domain. 
In Figure 1 1, we see the graphical construction of the coirelation-output for k = 0, N and 2N- 
5 1 . The convolution-results of k = 0 to k =■ N-l don't play any role, as the window after the 
IFFT in the PBFDAF discards these points. For the points N to 2N-1 in 

we see that convolution of tail 1 is largely concerned with input block B2, When we assume 
that this tail is small, this convolution can be approximated by the convolution of input block 
B2 with this tail in the next partition. This is achieved by removing tail 1 in the right side of 

10 the adaptive weights and placing this tail inside the first N points of the adaptive weight of 
the next partition, depicted in Figure 12. We see that for small tails, the convolution output 
concerning the original tail 1 in this partition is approximately equal to the convolution 
output of this tail in the next partition. 

On the other hand, we see that in Figure 12 there is also a right tail indicated 

15 by the number 2. For the points N to 2N-1 of the convolution output, we see that the 

convolution of tail 2 is also largely concerned with block B2. For convolution outputs close 
to 2N-1, there is only a small wrap-around with block B3. Again, for small tails, this 
convolution of tail 2 with input block B2 can also be obtained (approximately) in the 
previous partition. 

20 The destination of this tail 2 is depicted in Figure 1 1 . When the adaptive weights are shown 
with modified principal intervals as in Figure 9, we can visualize the compensation in the 
partition i as shown in Figure 1 3 . We note that this works well for the middle K-2 partitions, 
but there are problems for the first and the last partition. This topic is discussed later. 

Figure 14 shows a block diagram of an embodiment of a generator means for 

25 generating forward and backward compensation signals and Figure 15 shows some signal 
diagrams by means of which the operation of the generator means of Figure 14 is explained. 
For every weight vector constrained with a rectangular constraint, the forward and backward 
compensation signals 
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are generated. For reasons of transparency, all signals are constructed in time-domain, If we 
consider the unconstrained weight-vector as consisting of a linear convolution part and tailed 
circular wrap-around artifacts, the signals 1 to 6 of Figure 14 can be visualized as depicted in 
Figure IS. A PBFDAF using this compensation has an improved convergence curve which - 
5 does not have the large degradations in the NMSE on the moments that the exact rectangular 
constraint is performed. The positive effect of the compensation is clearly visible on the 
moments rectangular constraints are applied. The known PBFDAF shows large degradation 
in NMSE on these moments, whereas the compensated version shows an improvement. This 
improvement is due to the fact that the circular wrap-around artifacts are constrained to zero, 

10 while the useful information of the approximately linear correlation results is saved. 

Furthermore, the compensation of the circular-errors from one partition into the neighbor . 
partitions has a positive effect on the convergence time-constant or CTC. 

We compensate N/2 tune-domain points out of the circular wrap-around 
artifacts into the next and N/2 time-domain points into the previous partition. Only a part of 

15 these N/2 points are approximately linear (i.e. the tails in the figures) and thus useful, while 
the points between the tails are polluted with artifacts that have approximately 50% wrap- 
around, and can be assumed to be noisy and useless. When we also compensate die noisy 
artifacts as depicted in Figure 14, the steady-state NMSE is slightly worsened. A solution for 
this problem is nof to compensate the noisy part and only compensate the two tails. For this 

20 purpose, the compensated alternative constraining method in combination with the use of 
sinusoid windows as an approximation of the rectangular constraint may be used. The 
sinusoid window constrain the noisy parts in the update, while two tails remain for the 
compensation. 

Using an unconstrained PBFDAF, the convergence behavior is impaired and 
25 the region of convergence is smaller than with a fully constrained PBFDAF. When applying 
a sinusoid constraint approximation g 2 " in combination with the alternative constraining 
method, unproved results are obtained, compared with the unconstrained PBFDAF. The 



sinusoid approximation constraint g 21 * is defined as: 

When applying sintstsoid constraints, the noisy parts in the neighborhood of 3N/2 are reduced. 
30 On the other hand, the tails are not effectively reduced by the approximation constraint, but 
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22 21.08.2000 
can be successfully compensated into the neighbor partitions with the method described in 
the previous section. The application of the rectangular windows of length N/2 in Figure 14 
requires two extra FFTs. This complexity can be reduced by replacing the rectangular 
windows by a raised cosine and an raised inverse cosine window, as indicated in Figure 1 6, 
5 In Figure 1 4 the selection means comprise window block 104 and/or window block 106. 

Figure 1 6 shows a block diagram of another embodiment of a generator means 
for generating forward and backward compensation signals and Figure 17 shows some signal 
diagrams by means of which the operation of the generator means as shown in Figure 16 is 
explained. Again all signals in are constructed in time-domain for transparency reasons. The 
10 raised cosine windows can be easily implemented in frequency-domain without the use of 
extra FFTs. The sf^nals 1 to 6 are the same as in Figures 14 and 15 except that the noisy part 
is largely reduced by the sinusoid constraint approximations. The window for the backward 
compensation in the range i < 2N can be calculated by multiplying the sinusoid 
constraint window and the raised inverse cosine window. The window for the forward 
1 5 compensation can be calculated by multiplying the sinusoid constraint window qtiH the raised 
cosine window. Figure 18 shows the window for the backward compensation. Figure 19 
shows the window for the forward compensation. The selection means comprise here the 
combination of the sinusoid window and the raised cosine window (block 104)and/or the 
combination of the sinusoid window and the raised inverse cosine window (block 106). 
20 In figure 20 we can see that the compensation leads to connections between 

the update blocks 34 of adjacent partitions. For the first partition, there is no compensation to 
the previous partition. For the last partition we omit the compensation to the next (non- 
existent) partition. Omitting this forward compensation for the last partition leads to removal 
of circular wrap-around artifacts which could be useful for modeling a higher order impulse 
25 response. The update blocks 34 in Figure 20 exploit the compensation method as also 
explained in the previous section. The structure of these update blocks 34 is depicted in 
Figure 21. 

The sinusoid constraint 140 is placed before the coefficient update 142, while 
the rectangular constraint 150 is placed after the coefficient update 142, 144. The 
30 compensation signals are indicated with 



F^inted:OVO3-20d3 

- — — — _ — — — — — *£_-: — ^MJjt 



faWiM^ BBa 16: 13 " PHIlAcIP l li g- ' A~ «R.Z54 .-fffi^ 

PHNL000'475EPP 

» • 

23 21.08.2000 

AH {forward) and C^[x AT] (baclowird). 
When no rectangular constraint is applied the compensation signals are omitted (i.e. state 0 of 
the switches 146, 148, 152, 154). A linear way of rectangular constraining is applied, where 
the switches 146, 148, 1 52, 1 54 are described by the following formula: 



. if* mod = / 

0 otherwise * (5J> 



5 The use of an approximation constraint when also a rectangular constraint is 

applied seems to be needless. However, the compensation after applying a rectangular 
constraint also compensates the extra noise in the wrap-around artifacts (due to the omitting 
of the sinusoid constraint). Hence this leads to degradation in convergence behavior. 
Therefore it is preferred to also apply the approximation constraint when applying of a 

1 0 rectangular constraint with the compensation (see Figure 21). 

When implementing the rectangular constraint together with the 
compensation, we need additional computational complexity for the compensation part. This 
is because of the need of splitting up the two tails for the compensation in the next and 
previous partition and doing a time-domain shift of the two individual tails. A direct 

1 5 implementation as in Figure 14 is possible, but needs an abundant amount of extra 

computational complexity as all operations are performed in time-domain. A more efficient 
implementation can be obtained when generating the signals in frequency-domain. 
The shift in the time-domain is a N point shift of a 2N point signal 

Hi ♦ 

resulting in the 2h^point signal 
20 In frequency-domain these signals are represented by 
and 

respectively. This shift can be easily implemented in the frequency-domain via lemma 5.3.1 
which uses the shift-theorem. 
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Lemma 5.3.1 p 65? ^ea: 

<fl™\ * / fork*™* /<4X 

The generation of the time-domain signals 2, 3 and 4 of Figure 14 can be easily calculated in 
frequency-domain. Figure 22 shows this efficient implementation. We note that 

The ejementwise multiplication of a frequency domain vector by this vector implements the 
time-domain shift of N points. As the compensation signals need to be added to the neighbor 
5 partitions, these elementwise multiplications by -1 and 1 can be omitted if we use both 
additions and subtractions. The raised cosine window is now applied in frequency-domain. 
The other (raised inverse cosine window) is calculated by simple subtraction operations and 
saves computational complexity. It is also noted that the multiplications by 1/2 and j/4 can be 
implemented as simple shift operations, 

1 0 We already mentioned that it is not possible to perform the forward 

compensation in the last partition. Usually the circular wrap-around artifacts of this left tail 
are nulled by the rectangular constraint. As the wrap-around artifacts of this tail can be used 
for a higher order (>KN) adaptive filter, we could leave this tail unconstrained. The 
implementation of the compensation in the last partition is depicted in Figure 23. We note 

1 S that the computational complexity of this implementation is slightly lower compared to the 
implementation of Figure 22. 

In the embodiments of Figures 22 and 23, the selection means comprise the 
IFFT 160, zero appender 162, FFT 164 and the convolution block 168. 

We already mentioned that the first partition in the PBFDAF does not retrieve 

20 a forward compensation signal from its predecessor partition, and the first partition does not 
have to generate a backward compensation signal. Hence the implementation of the update- 
block in the first partition can also be implemented with reduced computational complexity. 

Good constraint approximations for the Partitioned BFDAF algorithm do not 
result in a near-optimal convergence behavior compared with the fully constrained PBFDAF, 

25 This is caused by unconstrained circular wrap-around artifacts, which are accumulated in the 
weight-update and cause problems in the convolution part of the PBFDAF as the different 
partitions are coupled. When combining the constraint approximations with the Alternative 
constraining mechanism, we get an improvement, but still achieve non-optimal convergence 
behavior. A specific part of the accumulated wrap-around artifacts not removed by the 
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25 21.08.2000 
constraint approximations are removed once in a while with the Alternative constraining 
mechanism. This specific part are approximately results of a linear correlation, that can be 
successfully reused (compensated) in the neighbor partitions. When we combine this method 
of compensation with the simplest form of constraint approximations (i.e, sinusoid 

- - - -5 constraint), we get a near-optimal convergence behavior compai^ to the fully^constrained _ 1 . _'■ 

PBFDAF, while obtaining an enormous reduction in computational 
complexity. For four and more partitions, the number of operations of the PBFDAF 
is reduced by a factor two. Simulations show that the performance of this constraining 
mechanism is also very good for highly correlated, non-stationary input signals. 

10 The PBFDAF according to the invention may be implemented in hardware, in 

software or in a combination of both. Preferably, a digital signal processor is used to 
implement the PBFDAF according to the invention. 

The scope of the invention is not limited to the embodiments explicitly 
disclosed. The invention is embodied in each new characteristic and each combination of 

1 5 characteristics. Any reference signs do not limit the scope of the claims. The word 

"comprising" does not exclude the presence of other elements or steps than those listed in a 
claim. Use of the word "a" or "an" preceding an element does not exclude the presence of a 
plurality of such elements. 
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Symbols 



a Mcsm value far constraint approximations, 

q Adapts lion constant. 

^/ Update block for partition;. 

A[ ■ } Numbor of addition*. 

Af J Backward compensation signal from partition /. 
C/Jk Afl Forward compensation signal from partition z. 

I>|> Rotation matrix. 

©{*} Number of divisions, 

^[*] Desired signal {echo in AEC). 

Estimate of the desired signal. 

5J*] Desired signal corrupted with 5 [4]. 

J* Fourier transform matrix. 

G 2r * Generalised! constraint window. 

Sl&l Stochastic gradient vector. 

Impulse response vector. 

1* identity matrix. 

j Imaginary imtt, 

J v Mirror matrix 

fC Number of partitions in the PBFDAF algorithm. 

K, Threshold partition for Hybrid constrained PBFDAF, 

k Block time index. 

m Multiplicative factor for high-slope constraint approximation. ; 

M[ * } Number of real multiplications* 1 

© Memoty occupation. { 

N Block fiength, \ 

g, O All zero vector, matrix, 

P Alternating constraint period. 

P x EA] Normalisation (power) vector. ^ 

Estimate of P^. 

* Computational complexity measure, 

r\k\ Residual signal. 

&[k] Near-end signal 
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FFT 

Invars* FFT 



Half Overlap 



Circular shift 



SerlaS-Parafllel converter 



Parallel/Serial converter 



l-to-2 switch 



2-to-l switch 



Append acros on the right 
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Nailing rlgbt part 



Dispose elements 



Raised Cosine window 
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Rectangular constraint 



Higfr-<>lope window 



2AT.<j 
0 

2N-a 



L 2 "l«*ri = ^(Jf^^AH ® G 2 *). where 



J for/ = 0 
22V- ly for/ =j 1 , 2iV — 1 
0 elsewhere 



£>"J = z* {X 2 * © where 



forimO 

0 "* fwO < / < even 

~ 2,V2t forO <: i < 2N, i odd 



for * =0 

for 0 < / t 2AT, / even 



IT- ^w- ^ - w-( 2a T i=i >^ j for 0 £ / < i odd 



Sinusoid! window 



Y 2 »lKtf$ a ^(^[eJV] where 

(2#-| for/ = 0 

2iV-i.2A'/ for/ = 2^-1 
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CLAIMS; 



L A partitioned block frequency domain adaptive filter for filteriiigah input 

signal in dependence on a control signal, the adaptive filter comprising a plurality of parallel 
arranged filter partitions, each filter partition being arranged for modeling a part of an 
impulse response of the adaptive filter, each filter partition having update means for updating 
5 filter coefficients of that filter partition by circular convoluting a signal representative of the 
input signal and a signal representative of the control signal, the update means comprising 
constraint means for intermittently constraining the filter coefficients by eliminating circular 
wrap-around artifacts of the circular convolution, characterized in that the update means 
further comprise selection means for selecting and removing at least part of the circular wrap- 
1 0 around artifacts, the selection means comprising an approximation of a rectangular constraint 
window. 

2. The partitioned block frequency domain adaptive filter according to Claim I, 
characterized in that the selection means are arranged for selecting the at least part of the 

1 5 circular wrap-around artifacts when the filter coefficients are not constrained by the 
constraint means. 

3. The partitioned block frequency domain adaptive filter according to Claim 1 
or 2, characterized in that time domain values of the approximation are larger than or equal to 

20 zero. 

4. The partitioned block frequency domain adaptive filter according to any one of 
the Claims I to 3, characterized in that the approximation in time domain has substantially 
high slopes near the positions which correspond to the positions of the transitions in a 

25 rectangular constraint window. 

5. The partitioned block frequency domain adaptive filter according to any one of 
Claims 1 to 4, characterized in that frequency domain values of the approximation each 
comprise a real value and complex conjugate imaginary values, whereby at least part of the 
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imaginary value&Jorm a row of numbers, the numbers being obtainable from one another by 
multiplication- 



5 



6. The partitioned block frequency domain adaptive filter according to Claim 5, 

characterized in that the frequency domain values of the approximation are defined as: 





^n-(^) - m-C^^)^ J for 0 < / < 2tf, i odd 



with i being an index number, m being a multiplication factor, a being a mean value. 



7. The partitioned block frequency domain adaptive filter according to Claim 6, 

characterized in that m is substantially equal to 2. 1 66. 



10 



8. 



An«acoustic echo canceller comprising a partitioned block frequency domain 



adaptive filter for filtering an input signal in dependence on a control signal, the adaptive 
filter comprising a plurality of parallel arranged filter partitions, each filter partition being 
arranged for modeling a part of an impulse response of the adaptive filter, each filter partition 

1 5 having update means for updating filter coefficients of that filter partition by circular 
convoluting a signal representative of the input signal and a signal representative of the 
control signal, the update means comprising constraint means for intermittently constraining 
the filter coefficients by eliminating circular wrap-around artifacts of the circular 
convolution, characterized in that the update means further comprise selection means for 

20 selecting and removing at least part of the circular wrap-around artifacts, the selection means 
comprising an approximation of a rectangular constraint window. 

9, The acoustic echo canceller according to Claim 8, characterized in that the 
selection means are arranged for selecting the at least part of the circular wrap-around 

25 artifacts when the filter coefficients are not constrained by the constraint means. 

10. The acoustic echo canceller according to Claim 8 or 9, characterized in that 
time domain values of the approximation are larger than or equal to 2ero. 

30 11. The acoustic echo canceller according to any one of the Claims 8 to 10, 

characterized in that the approximation in time domain has substantially high slopes near the 
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positions which correspond to the positions of the transitions in a rectangular constraint 
window. « 

12. The acoustic echo canceller according to any one of Claims 8 to 1 1, 
~ 5 characterized in that frequency domain values of the approximation each comprise a real 

value and complex conjugate imaginary values, whereby at least part of the imaginary values 
form a row of numbers, the numbers being obtainable from one another by multiplication. 

13, The acoustic echo canceller according to Claim 12, characterized in that the 
10 frequency domain values of the approximation are defined as; 



15 



2N>a for isO 

0 for 0 <. i < 2N, i even 

^ _ m-i^iF 1 )^ j for 0 < / < i odd 



"3" 

with i being an index number, m being a multiplication factor, a being a mean value, 

14. The acoustic echo canceller according to Claim 13, characterized in that m is 

substantially equal to 2.166. 



15. A method of adaptively filtering an input signal in dependence on a control 

signal, the method comprising the steps of: 

- partitioning the input signal into partitions, 

- for each partition updating filter coefficients by circular convoluting a signal representative 
20 of the input signal and a signal representative of the control signal, 

- for each partition intermittently constraining the filter coefficients by eliminating circular 
wrap-around artifacts of the circular convolution, 

characterized in that the method further comprises the step of selecting and removing at least 
part of the circular wrap-around artifacts by means of an approximation of a rectangular 
25 constraint window. 

The method according to Claim 15, characterized in that the selection of the at 
least part of the circular wrap-around artifacts is performed only when the filter coefficients 
are not constrained by the constraint means. 

30 
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17. The method according to Claim 1 5 or 1 6, characterized in that time domain 
values of the approximation are larger than or equal to zero. 

18. The method according to any one of the Claims IS to 17, characterized in that 
5 the approximation in time domain has substantially high slopes near the positions which 

correspond to depositions of the transitions in a rectangular constraint window. 

1 9. The method according to any one of Claims 1 5 to 1 8, characterized in that 
frequency domain values of the approximation each comprise a real value and/or an 

1 0 imaginary value and/or a conjugate imaginary value, whereby at least part of the imaginary 
values form a row of numbers, the numbers being obtainable from one another by 
multiplication. 



20. 



The method according to Claim 1 9, characterized in that the frequency domain 



2H a for i ^0 

ftiZNx _ 0 for 0 < ; < 2 A', i even 

-2^aI TC -l a ) - m -\— 5— ) \j f or 0 < i < 2N, i odd 

1 S values of the approximation are defined as: 

with i being an index number, m being a multiplication factor, a being a mean value. 



21 . The partitioned block frequency domain adaptive filter according to Claim 20, 

characterized in that m is substantially equal to 2. 1 66 . 
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ABSTRACT: 



_ ~ t^e partitioned bl^ 

invention comprises a plurality of parallel arranged filter partitions. Each filter partition 
models a part of an impulse response of the adaptive filter and has update means for updating 
filter coefficients*of that filter partition by means of a circular convolution. The update means 
5 intermittently constrain these filter coefficients by eliminating circular wrap-around artifacts 
of the circular convolution. The update means comprise selection means for selecting and 
removing at least part of the circular wrap-around artifacts. The selection means can be 
implemented with a relatively low computational complexity by means of an approximation 
of a rectangular constraint window. 

10 

Fig. 7 



Qni^ancLen) 21/ 8/0O 16:21 ; »31 4Q 27434B 9 -> EPO/EPA/OEB Rijswijk; Paglna 4Q 

21-O8-2OOo] 2000 16:20 ^^hilips cip r[EP0020291 6.3,39 ".' " NR. 254 -f[sPEC 



V 



Update 



A 



_ V 



■12. 




elk] 
(Echo) 



(Signal) 



Fig.l 



50% 

















«S*1 


_ 







Fig. 2 




3M 



in. 

L 



Fig. 3 



Printed:d1-63^26)0i 



6nt r yAnriRn; 21 / S/OO 16:21 ; 

PHIL] 



r2l5O8-2OOOi2000 16 : 21 



-+•31 _4jp_27^3^JB.9_-_>_ EPO/ EPA/ OEB R ± j swlj k ; Pagina 

NR. 254 



kCIP 



C: 



G*adfan*-conatciixit 
[ N 




5^ 



Fig. 4- 



^* Gxadisar*con&tiaint f 
*fV*1 . \^ i ' x 7 ~i "~~ i~ — " T ! . r-A Jg'l**} 

\ T *\ 




3f 



Fig. 5 




Fig. 6 




36? 




Oft 



Print|i:ll!03-200T 



_Orvt«_angen5 21/ 8/00 18:22; +ai tu -«g/«»o«*o » -- =■ — J — 

21-O8-2OOaf 0g0 16=21 ^Philips cip N' E poo20291 6.3; 3 A -nr.254 -fj SpEC 




Fig. 



Cataahr tUauff Circular 

Pitt pOB% 



0 AT 




Kg.**"-!? 



: i 



;2} 



1. 



C uts- 1 ^ 



3 ^ 



: > 



To.. 



:f 



64 



Pririted:01 .03^20C>1| 



38 



i:2i#8!ilPi2000 16=21 



PHIU 



ip mm%®2®m$b9 



CIP 



OEB Rij swiJK; Pagina -?r 



MR. 254 — I 



: f 



1 




» 4 3 

■ J <3> 



1 

5 



.•V 



RaJrad cevIim w&**r*- 



2JV 



mmmmmmm 



P39J 



<S r».*-_V« r> o *> n ^ 21/ B/OO 16:22; 

21 -08-2000 ]2000 16:21 ^^HILIPS CIP ^EJf^|^^|p|@^ 9 



+ 31 ^O^T^^eg^^ EPO/EPA/OE8 Rijswijki PagAne 4 



R.B£f ad Envtirfia ctsgM window 




Fig. yf 17 



r h 








ll 



3© 



Aq Si i iri s \i ii , = aaa^ 



3 

Sn 




i r 

6ntVangen* 21 / 8/00 16:22; 

pj^||Pig|28B9 16:21 PHIL 



+31 40 8743489 -> EPO/EPA/OEB Rljswi]k; Paglna 47 



^cip r,,E.| 



te^0Of,©29;1t6i'3*39 



NR. 254 ~ H 




CyfV-VJ 



lD 



T 



-At 




nc 

It . "«t 



=0 





=Ir v 

> 



Fig. £3"^ 



2W, 



T 
It*. 



t 

11.V 



I 



1 i«t? 1 L 1 



Nnted:01-03-2Q(ft 



1 ^ 



CERTIFIED COPY OF 
PRIORITY DOCUMENT 



